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Abstract. We present a high-resolution spectroscopic study of the envelope 
of the young Class protostar IRAM 04191+1522 in Taurus. N 2 H+(l-0) 
observations with the Plateau de Bure Interferometer and the IRAM 30m 
telescope demonstrate that the molecular ion N2H 4 " disappears from the gas 
phase in the inner envelope (r < 1600 AU, n H2 > 4 — 7 x 10 5 cm -3 ). This may 
result from N2 depletion on polar ice mantles and enhanced grain chemistry. 



1 The Class Protostar IRAM 04191 

With an age t ~ 1 — 3 x 10 4 yr measured from the beginning of the accretion 
phase, IRAM 04191+1522 is the youngest Class protostar known so far 
in the Taurus molecular cloud (d = 140 pc). It has a prominent envelope 
(1.5 M Q ) and an associated powerful bipolar outflow [3]. Belloche et al. [4] 
showed that the envelope is undergoing both extended infall motions and 
fast, differential rotation. They proposed that the rapidly rotating inner part 
of the envelope (r < 3500 AU) corresponds to a magnetically supercritical 
core decoupling from an environment still supported by magnetic fields and 
strongly affected by magnetic braking. 

We recently carried out new observations in the N2H + (l-0) line with the 
Plateau de Bure Interferometer (PdBI) to probe the inner structure of the 
IRAM 04191 envelope. Here, we present and discuss the results of these high- 
resolution (~ 5") observations of the small-scale distribution of N2H" 1 ". 



2 Small-Scale N 2 H+(l-0) Emission 

Fig. la shows the N2H + (l-0) integrated intensity map obtained with PdBI, 
which filters out the extended emission observed with the IRAM 30m 
telescope. The (naturally-weighted) synthesized half-power beamwidth is 
5.8" x 4". Two emission peaks are seen approximately along the direction 
perpendicular to the outflow axis, on either side of a gap very close to the 
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Fig. 1. (a) N 2 H+(l-0) integrated intensity map of the IRAM 04191 envelope ob- 
tained with PdBI. The contour step is 0.1 Jy beam -1 km s -1 (i.e. 3 times the rms 
noise) . The big circle represents the primary beam. The emission gap very close to 
the central position (black cross) suggests a strong decrease of the N2H + abundance 
in the inner part of the envelope, (b) Map of the LSR velocity (squares) as derived 
from Gaussian hfs fits to the N2H + (l-0) spectra. The map of integrated intensity 
is overlaid as contours. We measure a mean velocity gradient of ~ 26 km s _1 pc _1 
across the envelope (black arrow), whose direction at PA ~ 107° is consistent with 
the rotation seen on larger scale. 



central position. This suggests a strong decrease of the N 2 H + abundance 
toward the center. Besides, the two emission peaks have very different cen- 
troid velocities (see Fig. lb): we measure a strong velocity gradient of ~ 26 
km s _1 pc" 1 along the axis with PA = 107°, i.e. approximately the same 
axis along which a rotational velocity gradient is detected on larger scale [4] . 
The gradient measured with PdBI is 8 times larger than the mean velocity 
gradient observed with the 30m telescope toward the core, which confirms 
the presence of fast, differential rotation in the IRAM 04191 envelope. 



3 N 2 H+ Disappearance from the Gas Phase 

Since the interferometer filters out the extended emission, we added short 
spacings obtained with the 30m telescope [4] to the PdBI data. As shown 
in the N2H + (l-0) integrated intensity map (see Fig. 2a), the emission is re- 
solved in a ringlike structure with two local peaks, surrounding a gap close 
to the center. Gaussian hfs fits to the 7-component N2H + (l-0) spectra yield 
a mean total opacity Ttot = 10 ± 4 in the central part (9 < 10"), correspond- 
ing to r = 2.6 ± 1.0 for the strongest component and r = 1.1 ± 0.4 for the 
isolated component. However, the opacity map does not show any enhance- 
ment toward the center and the gap close to the center does not have the 
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Fig. 2. (a) N2H + (l-0) integrated intensity map of the IRAM 04191 envelope com- 
bining the PdBI data and the short spacings obtained at the 30m telescope. The 
contour step is 0.16 Jy beam -1 km s -1 (i.e. 7 times the rms noise). The white circle 
represents the PdBI primary beam, (b) Radial profile of the N2H + (l-0) integrated 
intensity corrected for the primary beam attenuation. The dots with error bars rep- 
resent the average intensity in concentric shells limited to the positions shown in the 
insert. The solid curve indicates the H 2 column density profile [10, 4] scaled to the 
9 th shell for comparison. The dashed curves show the effect of a hole of N2H + (l-0) 
emission with radius r = 1350 and 1850 AU, respectively. They suggest a strong 
decrease of the N2H + abundance in the inner part of the envelope (r < 1600 AU). 

highest opacity. Besides, the intensity maps integrated over each AF\ group 
of N 2 H + (l-0) hyperfine components separately show the same pattern. We 
are thus confident that the observed emission gap docs not result from an 
opacity effect. 

Assuming a constant excitation temperature for N2H + (l-0), we compare 
in Fig. 2b the circularly-averaged profile of the N 2 H + (l-0) integrated in- 
tensity, corrected for PdBI primary beam attenuation, with the H2 column 
density profile derived from the dust emission [10, 4]. The N 2 H + (l-0) profile 
departs from the H 2 column density profile for 9 < 13". A good model fit is 
obtained when one assumes a hole of N 2 H + (l-0) emission in the central part 
of the envelope (r < 1350-1850 AU, see dashed curves). This crude compari- 
son thus strongly suggests that the N 2 H + ion disappears from the gas phase 
above a density of n H2 ~ 4 - 7 x 10 5 cm -3 in the IRAM 04191 envelope. 

4 Discussion 

Very little evidence for N 2 H+ abundance drop has been reported so far in 
dense cores (see however the case of B68 [6]). The N 2 H + (l-0) emission usu- 
ally follows the dust emission very well in prestellar cores and has long been 
thought to be one of the best molecular tracers of dense gas [12, 7]. From 
a theoretical point of view, all recent chemical models predict a centrally- 
peaked N 2 H + column density profile during the prestellar collapse phase 
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[5, 1, 2, 9, 11]. In particular, the models calculated by Aikawa et al. [2] in- 
cluding gas-phase reactions, gas-dust interactions and diffusive grain-surface 
reactions predict an increase of the N 2 H+ abundance during the core collapse 
up to nn 2 ~ 1 x 10 7 cm~ 3 or more. Since only a few 10 4 yr have elapsed and 
the density at r ~ 1600 AU cannot have changed much between point mass 
formation (the last stage of this model) and the present stage of IRAM 04191, 
the drop observed in the N 2 H + abundance above nn 2 ~4-7x 10 5 cm~ 3 is 
strongly at variance with current theoretical predictions. 

Our result implies either a stronger dependence of the N 2 H + abundance 
on density (i.e. an earlier drop in abundance during the prestellar phase) or 
an abundance drop on a timescale as short as a few 10 4 yr. Since N 2 is the 
precursor of N 2 H + , the drop in N 2 H + abundance may result from stronger 
N 2 depiction onto grain surfaces. This could happen if N 2 were depleted on 
H 2 ice mantles characterized by a higher binding energy (see [6], Model 
C of [1] and Fig. 6 of [11]) and/or if the chemistry transforming N 2 in less 
volatile species on grain surfaces were enhanced. It is also possible that N 2 H+ 
destroyers such as CO have returned to the gas phase in the inner envelope, 
following heating by the central protostar or shock by the outflow. However, 
the dust temperature is below 10 K at r ~ 1000 AU and our C 18 0(2-l) 
observations [4] at a resolution of 800 AU in radius show depiction toward 
the envelope center by a factor of 3 in column density when compared to the 

I. 3mm continuum emission. This latter explanation seems thus less likely. 

A better tracer than N 2 H + will thus have to be found to probe the veloc- 
ity structure of cold protostcllar envelopes at high resolution with the next 
generation of (sub)millimeter instruments such as ALMA. H 2 D+ may be such 
a good tracer [6, 8]. 
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